Dark energy, matter creation and curvature 
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The most studied way to explain the current accelerated expansion of the universe is to assume 
the existence of dark energy; a new component that fill the universe, does not clumps, currently 
dominates the evolution, and has a negative pressure. In this work I study an alternative model 
proposed by Lima et al. [Sj, which does not need an exotic equation of state, but assumes instead 
the existence of gravitational particle creation. Because this model fits the supernova observations 
as well as the ACDM model, I perform in this work a thorough study of this model considering 
an explicit spatial curvature. I found that in this scenario we can alleviate the cosmic coincidence 
problem, basically showing that these two components, dark matter and dark energy, are of the 
same nature, but they act at different scales. I also shown the inadequacy of some particle creation 
models, and also I study a previously propose new model that overcome these difficulties. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

Currently the observational evidence coming from su- 
pernovae studies [l|, cosmic background radiation fluc- 
tuations 0, and baryon acoustic oscillations 0, set a 
strong case for a cosmological (concordance) universe 
model composed by nearly 70 percent of a mysterious 
component called dark energy, responsible for the current 
accelerated expansion, nearly 25 percent of dark matter, 
which populates the galaxy halos and a small percentage 
(around 4) is composed by baryonic matter. The na- 
ture of these two dark component remains so far obscure 
In the case of dark matter, we have a set of candi- 
dates to be probed with observations and detections in 
particle accelerators, however the case of dark energy is 
more elusive. This component not only has to fill the 
universe at the largest scale homogeneously, and also has 
the appropriate order of magnitude to be comparable to 
dark matter, but also has to have a negative pressure 
equation of state, something that was assumed first to 
explain inflation in the early universe, inspired by high 
energy theories of particle physics, but which seems to be 
awkward to appeal for, at these low energy scales. The 
alternative way to account for the cosmic acceleration, 
in the framework of the standard model, is to consider a 
modification of general relativity at large scales 

Some years ago, Prigogine and co-workers @ presented 
a very interesting cosmological model where matter cre- 
ation takes place without spoiling adiabatic expansion. 
This is possible by adding a work term due to the change 
of the particle number density n. This work actually 
suggested for the first time, a way to incorporate the 
particle creation process in the context of cosmology in 
a self consistent way. In fact, the original claim made by 
Zcldovich [7J, that gravitational particle production can 
be described phenomenologically by a negative pressure, 
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is realize here in a beautiful way. A covariant formulation 
of the model was presented for the first time in 

Originally, these considerations were applied in the 
context of; the steady-state cosmological model, the 
warm inflationary scenario and within the standard infla- 
tionary scenario, during the reheating phase. Based on 
this workvthe studies of cosmological models with matter 
creation |9| were initiated, and rapidly reco gniz ed to be 
potentially important to explain dark energy [lfj . In par- 
ticular, in Rcf.[ll] the authors established a model where 
dark energy can be mimicked by self-interactions of the 
dark matter substratum. Within the same framework, 
models of interacting dark energy and dark matter were 
proposed [liq . Actually, in these studies is possible to 
have consistently, a universe where matter creation pro- 
ceeds within an adiabatic evolution. More recently, Lima 
et al.,[H| have presented a study of a fiat cosmological 
model where a transition from decelerated to accelerated 
phase exist. They explicitly show that previous mod- 
els considered jlif . does not exhibit the transition, and 
study the observational constraints on the model parame- 
ters. Also in the authors modified slightly the model, 
adding explicitly a baryonic contribution, which enable 
them to have a transition from decelerated to accelerated 
expansion always. 

In this work I consider the non flat extension of this 
matter creation cosmological model. This is important 
because, even in the ACDM concordance model, there is 
no clear evidence that flk is zero, due to the well known 
degeneracy between w(z) and fl6[. Furthermore, any 
alternative model to dark energy then must consider a 
non flat assumption as a prior. This model enable us 
to explain the current acceleration of the universe ex- 
pansion through a fictitious pressure component coming 
from changes in the dark matter particle number, with- 
out any exotic contribution. As a bonus, this model en- 
able us to explain easily the cosmic coincidence prob- 
lem; it is not strange to have a similar contribution from 
these two (commonly differently regarded) components, 
because they are just two aspects of one and the same 
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component; the dark matter. 

In the next section I derive the equations of motion 
in the case of matter creation. Then, in section III I 
study the case for non flat models and the observational 
consequences in the models already knew. After this, 
I discuss a new particle creation model that resembles 
the ACDM model case, showing the transition from a 
decelerated to an accelerated expansion without enter in 
conflict with the behavior for large z. 

II. THE EFFECTIVE NEGATIVE PRESSURE 

Assuming that the particle number is not conserved 
Nit 7^ 0, it leads to a modification of the energy conser- 
vation equation. In fact, assuming Nf^ = nT where T 
is the particle creation rate, and using the FRW metric 
we obtain the generalization of the energy conservation 
equation, 

d(pV) + pdV - (h/n)d(nV) = 0, (I) 

where h = (p + p) is the enthalpy (per unit volume) , n 
is the number density and p is the energy density, and a 
new equation 

n + 3Hn = nT, (2) 

for the non conserved number density. The important 
thing to stress here is that the universe evolution con- 
tinues to be adiabatic, in the sense that the entropy per 
particle remains unchanged (s = 0). The extra contribu- 
tion can be interpreted in ((T|) as a non thermal pressure 
defined as 

This is the source that produces the acceleration of the 
universe expansion. Once the particle number increases 
with the volume, we obtain a negative pressure. 

Now, because we are considering matter creation, we 
have to impose the second law constraint 

dS = -d(nV) > 0, (4) 
n 

where s — S/V is the entropy density. From (fTJ) we find 
that the new set of Einstein equations are 

n k 8itG ,_. 
H+ V=— P (5) 

Tl 

P=-{P+P) (6) 
n 

and the previously derived relation ([2]). The set of equa- 
tions ((2j, (J5j> and ([6]) completely specified the system 
evolution. It is also useful to combine (JU) with <j5j> and 
([5]) to eliminate p and obtain 



The standard adiabatic evolution is easily recovered: set- 
ting r = implies that h/n — —3H, which leads to the 
usual conservation equation from (|6]). A class of de Sit- 
ter solution is obtained with h = p = and arbitrary 
pressure p. Moreover, there exist another class of solu- 
tions where Eq.(j6]) enable us to determine the pressure; 
for example if p = mn, with m constant, Eq.(j6]) implies 
p = 0, and furthermore if p — aT and n = bT 3 implies 
V = p/3- 

I have to stress here that the main result derived in this 
section means that we have a single contribution, which 
satisfy the non-relativistic matter equation of state, that 
describe both dark matter and dark energy simultane- 
ously, and in this way solves automatically the coinci- 
dence problem. This dark unification mechanism does 
not have the problem studied in [2(| , because in the cases 
studied in that paper, for example the Chaplygin gas 
model jU] , the sound velocity of the dark matter is not 
zero, leading to instabilities. That happens because the 
Chaplygin gas interpolates between the equation of states 
for dark matter and dark energy. This does not happens 
here, because there is just one equation of state; that of 
non-relativistic matter. 



III. MATTER CREATION IN A NON FLAT 
UNIVERSE 

In this section we study a number of models where 
matter creation produces cosmic acceleration, in a curved 
background. We also test the models using the most 
recent Supernovae data, the so called Union 2 set [22[ • In 
order to do that test we consider the comoving distance 
from the observer to redshift z is given by 

V-^fc Jo E i z ) 

where E(z) = H(z)/Hq. The SNIa data give the distance 
modulus p(z) which is related to the luminosity distance 
d L (z) = (l + z)r(z), through p, (z) = 51og 10 [d L (z)/Mpc] + 
25. We fit the SNIa with the cosmological model by min- 
imizing the x 2 value defined by 

2 _ [p{Zj) - p ohs {zi)] 2 
XSNIa — 2 , 2 ' W 

where p, bs is the corresponding observed one. 



A. P Model 

Let us assume first the particle creation rate evolving 

as 

T = 3#j8, (fO) 
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where (3 is a constant. From ([2]) we finds the solution for 
the density number 
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3(1-/3) 



(11) 



where no is a constant of integration. Assuming p = mn 
as for non-relativistic matter (that leads to que usual 
equation of state p — as we discussed that at the end 
of the last section) , we get for the Hubble function 

e 2 (z) = n m (i + z) 3 ' 1 -^) + (i - n m )(i + zf. (12) 

Clearly, this model does not describe properly a transi- 
tion from decelerated to an accelerated expansion phase. 
In fact, from ([5]) and © we obtain: 



(13) 



Depending on the value assumed for /3 we obtain a model 
that accelerate forever ((3 > 1/3) or decelerate forever 
(/3 < 1/3). A Bayesian analysis using the Union 2 data 
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FIG. 1: Using the Union 2 data set [22j we plot the confidence 
contours, at 68.27% and 95.45% for the two parameters of the 
model: /3 and Qk- 

set of SNIa 22] enable us to plot confidence contours for 
the two parameters of the model: /3 and f2fc. This is 
shown in FigJTJ Clearly, the data suggest a large nega- 
tive value for the curvature parameter, indicating that in 
order to fit the SNIa data, the function H (z) has to have 
a very distorted form. In any case, because this model 
does not present a transition from a decelerated to an 
accelerated phase, it can be ruled out immediately. 



on what is used to model dark energy, different redshift 
have been obtained between 0.5 to 1. In the context 
of a model of adiabatic matter creation, this topic was 
recently discussed in 13[ where the authors specialized 
in a flat cosmology, where a explicit transition can be 
achieved from a decelerated to an accelerated expansion 
using the following model 



r = 3 7 if - 



(14) 



In this section we generalize this work to non flat uni- 
verses. In what follows, a non relativistic matter equa- 
tion of state is assumed p = mn (p = 0). 

Introducing (fl~4"|) in ^ leads to the following solution 



( ^ ) , :i ""<" '"' 



(15) 



Is evident the meaning of the subscript zero. Because p = 
mn for non-relativistic matter and using we obtain 



R , r ( H 1 



(16) 



that clearly allow us to describe an accelera- 
tion/deceleration transition. Writing the right hand side 
in terms of H and using the Hubble equation once more, 
we can write 



2 1 R ) \ H 3 



(17) 



For k = this equation coincides with that in [13j, for 
which an explicit form of H (z) was found. In the general 
case, k ^ 0, it is not possible to integrate (fTTj) to obtain 
a closed analytic form. Instead, we use R/R = H to 
change the time derivatives for derivatives of the scale 
factor R, and using that Rq/R = 1 + z, change again 
the derivatives in terms of the redshift. Doing that, we 
obtain 



rIE 

(l + z)E^--E 2 = [n fc (i 



2E 2 



(18) 

where E(z) = H(z)/Hq. Given the set of parameters 
fife and 7, we integrate numerically (1181) to obtain the 
function E(z). A Bayesian analysis using the Union 2 
data set of SNIa [22j enable us to plot confidence contours 
for the two parameters of the model: 7 and fifc. This is 
shown in Fig[2] Again, as in the previous subsection, the 
SNIa data suggest a very large negative value for il^ . 



B. 7 model 

As is well known from observations of supernovae, a 
transition from a decelerated to an accelerated expansion 
occurs in the recent history of the universe. Depending 



C. Problems with these models 

A simple way to understand the large negative values 
obtained for in the fitting process, which is suggested 
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FIG. 2: Using the Union 2 data set [22] we plot the confidence 
contours, at 68.27% and 95.45% for the two parameters of the 
model: 7 and Q^, by integrating Eq. QlSp . 



by the SNIa data in these models, is by plotting the H (z) 
function in the range z £ (0, 4). In Fig. [3] we display such 
a plot. All the three models have two free parameters to 
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FIG. 3: The Hubble function in terms of the redshift for the 
/3 model (dashed line), the 7 model (dotted line) and the flat 
LCDM model, (fi m ,n A ) = (0.27,0.73) (continuos line). 

be fixed with the SNIa data. Although at low redshift 
(z < 0.2) the three curves coincide, their evolution with 
increasing z differs appreciably. Clearly the f3 model be- 
haves erroneously. Instead, the 7 model, although it has 
the correct tendency, it is easily differentiated from the 
LCDM at large redshift (z > 0.2). Because we are using 
the SNIa data set, the best fit values for the free param- 
eters in these models (shown in FigEJ) adjust principally 
the low- 2 zone, because from the 557 SNIa data only 19 
have z > 1. This is also clear once we try to perform 
a joint Bayesian analysis using SNIa plus the BAO and 
CMB constraints (see Table I). Using the f3 model, the 
procedure even does not find any appropriate configura- 
tion. Using the 7 model, it is possible to get a good fit 
{Xred = 0.994) using only SNIa+BAO, however once we 
add the CMB constraint, it worsens the fit (xred = 1-32). 



TABLE I: The \ 2 f° r the best fit values in models discussed 
in this section plus the curved ACDM using the Union2 data 
set in the case of a non flat universe model. 



Model 


SN 


SN+BAO 


SN+BAO+CMB 


P Model 


542.17 


554.99 




7 Model 


542.29 


552.66 


731.12 


ACDMk 


542.55 


542.65 


542.65 



It is therefore clear why in [15] the authors add an 
explicit contribution for baryons. Adding baryonic dark 
matter with = 0.042 in the previous models, enable 
us to obtain better best fit parameters for the models 
that, for instance modifying Fig. [3] into Fig. |4] 
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FIG. 4: The same as Fig. O but this time considering bary- 
onic dark matter with Qb = 0.042. See the text for a discus- 
sion. 

For example, in the case of the 7 model, once we add 
baryons the best fit parameters improve respect to the 
case without baryons; using only SNIa+BAO we get a 
Xred = 0.989, however once we add the CMB constraint, 
it worsens the fit again. 

IV. SEARCH FOR A NEW MODEL 

The previous models have been discussed in several 
papers, deriving their consequences in cosmological evo- 
lution and also testing their performance to fit the su- 
pernova data. The problem with the (3 model is that it 
can not describe properly a transition from a decelerated 
to an accelerated phase. In contrast, the 7 model can 
do this, however they can only fit well low-redshift ob- 
servational constraint (SNIa data with z < 0.2). As we 
mentioned in the last section, if we attempt to include the 
BAO and CMB constraint, the whole fit breaks down. 

The objective merit of this type of model is that we 
can mimic a cosmological constant without use of an ex- 
otic component, we just need to assume the existence of 
matter creation that takes place with a certain rate. In 
principle, we do not have any fundamental reason to use 
a specific form of the matter creation rate. The only fun- 
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damental feature we have to respect is the second law of 
thermodynamics (QJ. 

In [T3| we propose a model that was also discussed in 
Let us consider this model here as an example. This 
model is characterized by the following matter creation 
rate 



nT = 3aH, 



(19) 



which for a > satisfy the requirement, because H (z) is 
a positive definite function. Inserting in ([2]) we obtain 



R 3 

n(R) = (n - oO^f + a - 



(20) 



that resembles the combined contribution of a cosmo- 
logical constant and dust. The level of fine tuning here 
to obtain an accelerated expansion is relatively smaller 
than in the usual ACDM model, because in this case, 
both terms come from the same function, and we know 
that we can obtain an accelerated expansion phase after 
some period of time from ([3]). In fact, from Q we get 
p c = —ma, so replacing in (JT)) we find 



R 
R 



AnGm 



2a — (no — a) 



R 



(21) 



Note however that the equation of state has not been 
modified, p is still the non-relativistic contribution. We 
do not have to introduce any exotic component - with 
a negative pressure - to describe the current expansion 
acceleration. Given the current status of the dark matter 
and dark energy problem we can use this idea as a 
possible way to understand it. Replacing in the Hubble 
equation we obtain 

E{zf = Q a + (1 - fl a - n*)(l + zf + n fc (l + zf, (22) 

where il a = ma/p cr and p cr = 3Hq/8ttG, which is ef- 
fectively indistinguishable from a curved LCDM model. 



Because ft a — = 0.7, is clear that the observations 
implies a positive a for this new model. 

In this work I have studied a model which considering 
small changes in the total number of particles in our uni- 
verse, offer a possible way to understand the current ac- 
celerated expansion measurements, without using any ex- 
otic energy component. Also, this scenario alleviate the 
cosmic coincidence problem, basically showing that these 
two components, dark matter and dark energy, are of the 
same nature, but their act at different scales. This way of 
understand the SNIa observations, implies that cosmol- 
ogy have a new window to explore the universe consid- 
ering matter creation, encoded in the function dN/dV. 
I have discussed the two models already known in the 
literature, extending the analysis to nonflat universes, 
performing a bayesian analysis using both SNIa and also 
BAO and CMB. I have found that, both the /3 and 7 mod- 
els fail to fit the observations, even considering fi^ 7^ 0. 
Although adding an explicit baryonic dark matter contri- 
bution may ameliorate the fit using SNIa+BAO in the 7 
model, adding the CMB constraint worsens the situation. 
The statistical analysis was performed even in the case 
an explicit form for H(z) was lacking. I have also stud- 
ied a previously proposed new model, through a specific 
form of the matter creation rate, which is almost indis- 
tinguishable from a LCDM model, which enable us to fit 
all the observations, SNIa+BAO+CMB. Although sim- 
ilar to LCDM, the model can be distinguished through 
the structure formation process, which will be studied 
elsewhere. 
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